Spin dynamics in a hole-doped 5 = 1/2 Heisenberg antiferromagnet with a disordered 

ground state 



(N 
O 

o 

(N 



Wei Bao and J.L. Sarrao 
Los Alamos National Laboratory, Los Alamos, NM 87545 
(Dated: February 1, 2008) 

Only 3% hole doping by Li is sufficient to suppress the long-range antiferromagnetic order in 
La2Cu04. Spin dynamics in such a disordered state was investigated with measurements of the 
dynamic magnetic structure factor S(tJ, q), using cold neutron spectroscopy, for La2Cuo.94Lio.06O4. 
The S(u), q) is found to sharply peak at (7r,7r), and its dynamics to be relaxational. Confirming 
theoretical expectation for the quantum disordered 2D S — 1/2 Heisenberg antiferromagnet, the 
energy scale saturates at a finite value at low temperatures. Possible connection to the "pseudo spin 
gap" phenomenon observed in the NMR/NQR studies on underdoped cuprates is discussed. 
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Stimulated by the discovery of cuprate superconduc- 
tors, which derive from doping charge carriers into weakly 
coupled two-dimensional (2D) antiferromagnetic Cu02 
planes, there has been great interest in 2D 5=1/2 Heisen- 
berg antiferromagnet (HAF) with dominant nearest- 
neighbor exchange interactions on a square lattice. It 
is now generally accepted that there is long-range Neel 
order at T = for such a system ffl. Finite-temperature 
magnetic properties Q are in good agreement with theo- 
retical predictions ||. 

The situation more closely related to superconductiv- 
ity in cuprates, namely, 2D 5=1/2 HAF with doped 
holes, however, is much less understood. The parent com- 
pounds such as La2Cu04 are charge-transfer materials^ , 
therefore, the doped hole is a charge with 5 = 1/2 
located at an O site. However, the formation of the 
Zhang- Rice singlet || allows an effective description of 
the hole as a spinless charge at a Cu site. Strong sup- 
pression of the Neel order by holes can be accounted 
for by long-range topological disturbances accompany- 
ing holes In addition, hole motion is also disruptive 
to the Neel order Q. For hole-doped La2CuC>4 with- 
out the Neel order, two type of incommensurate anti- 
ferromagnetic correlations have been discovered experi- 
mentally when holes introduced by dopants such as Sr 
or Ba are mobile||. When holes are loosely bound in 
the case of Li-doped La2Cu04|p|, antiferromagnetic cor- 
relations remain commensurate [ 1 Of] - Microscopic under- 
standing of interaction between holes and spins on the 
square lattice remains a major challenge in condensed 
matter research |pr| [l2"f . 

An alternative approach to 2D 5 = 1/2 HAF with 
doped holes is based on quantum phase transition. Build- 
ing on the success of the quantum non-linear a model 
as an effective low temperature theory for 2D 5=1/2 
HAF|3|, [l3|, effect of doping is simulated by a frustra- 
tion parameter g, which at a critical value g c suppresses 
the long-range Neel order at T = 0|| O. The advan- 
tage of this approach is that aspects of spin dynamics 
can be predicted from general theoretical arguments for 
dynamic critical phenomena before microscopic theory is 



established. The magnetically disordered state in the 
g-T plane is divided into three physical regimes: the 
renormalized classic, quantum critical (QC) and quan- 
tum disordered (QD). For doped cuprates with disor- 
dered ground state (g > g c ), only the QC and QD regimes 
are relevant. Quantum critical theory 1 14 |l5[ predicts 



that the energy scale, T, is proportional to T in the QC 
regime at higher T. Inelastic neutron scattering studies 
on Ba or Sr doped La2Cu04 have confirmed this pre- 
diction over an extraordinarily wide energy (hu> < 90 
meV) and temperature (T < 500 K) range 0. At 
lower T in the QD regime, T is expected theoretically 
to saturate at a finite value, Tq. Here we report a di- 
rect observation of such a behavior in a cold neutron 
inelastic neutron scattering study on La2Cuo.94Lio.06O4. 
With an energy resolution of 0.1 meV in term of the full- 
width-at-half-maximum of incoherent scattering, we fo- 
cused on low energy spin dynamics, relevant in search for 
To. Possible connection of the QD spin dynamics to the 
so-called "pseudo spin gap" phenomenon observed in the 
NMR/NQR studies on hole-doped La2Cu04 is discussed. 

A single crystal of La2Cuo.94Lio.06O4, weighing 2.1 g, 
was grown in CuO flux, using isotopically enriched Li 
(98.4%) to reduce neutron absorption. The crystal has 
orthorhombic Cmca symmetry with lattice parameters 
a = 5.3511, 6 = 13.151 and c = 5.3861 at 295 K. Using 
this orthorhombic unit cell to label reciprocal q space, the 
(tt, 7r) point in the square lattice notation splits to (100) 
and (001) points. Neutron scattering signal is observed at 
(100) type but not at (001) type Bragg points, as in the 
stoichiometric antiferromagnet La2Cu04|l^|. Measure- 
ments of seven independent (100) type of Bragg peaks, 
using thermal neutrons to reach 6 A , confirm the mag- 
netic origin of these peaks. In the remainder of this study, 
we will focus on spin dynamics near Q=(100), using cold 
neutron triple-axis spectrometer SPINS at NIST with 
fixed Ef = 3.7 meV or E f = 5meV. The (002) reflection 
of pyrolytic graphite was used for both the monochroma- 
tor and analyzer. A cold Be or BcO filter was used to 
eliminate higher order neutrons. Horizontal Soller slits of 
80' were used before and after the sample. Temperature 
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FIG. 1: (a) Const. -So; scans with various energy transfers 
at 30 K. (b) Constant-ftoj = 0.6meV scan at various tem- 
peratures. The q=(100) here in the orthorhombic notation 
corresponds to the (it, n) point of the Cu02 square plane. 



of the sample was regulated by a pumped He cryostat. 

Fig. ^ shows some constant-energy (hco) scans around 
Q=(100) at various energies and temperatures, which 
roughly cover the energy and temperature range of this 
study. All the peaks at (w, if) are resolution-limited. 
This means that dynamic spin clusters with nearest- 
neighbor antiferromagnetic alignment in the Cu02 plane 
have grown to substantial size. From the half-width-at- 
half-maximum of the peaks, a lower bound of magnetic 
correlation length can be estimated, namely, £ > A2A be- 
low ^100 K. The correlation length is much longer than 
the mean distance, 15A, between Li dopants. Therefore, 
consistent with expectation from microscopic theory by 
Haas et al.Q, antiferromagnetic correlations in our sam- 
ple are not simply impurity limited. 

Dynamics of the antiferromagnetically correlated spin 
clusters in La2Cuo.94Lio.06O4 are probed by energy scans 
at Q=(100) at various temperatures (Fig. ||). Since mag- 
netic intensity is confined within resolution at (100), 
the energy scan at (1.39,0,0), which is T and Tiuo- 
independent (refer to triangles in Fig. ||), offers a good 
measure of background. Lowering temperature from 
150 K, magnetic intensity at lower energies increases 
at the expense of magnetic intensity at higher energies, 
as shown by filled and open diamonds. This is typical 
for paramagnetic fluctuations, as magnetic susceptibility 
increases while fluctuation energy decreases with lower- 
ing temperature ]l9| . Upon further lowering temperature, 
however, magnetic intensity below ~ 1.5meV is progres- 
sively suppressed. At 1.5 K, a peak at a finite energy can 
be clearly discerned. We would like to emphasize here 
that there is no real gap in the spin excitation spectra 
in Fig. ||. As will be shown later, spectra in Fig. || can 
be well described by relaxational spin dynamics. The- 
oretically, both gapped and gapless spin dynamics have 
been found in different toy models in the QD regime G3. 
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FIG. 2: Const.-Q=(100) scans from 1.5 to 150 K. The solid 
curves at finite Tiuj and their dot-line continuation near lo = 
are theoretical fit to Eqs. (Q) and (^|) with parameters T 
and \Q shown in Fig. |j| Triangles represent a measured flat 
background at (1.39,0,0). 



The spectra in Fig. || also differ from inelastic neutron 
scattering spectra from conventional, heavy fermion, and 
cuprate superconductors, which do show spin or pseudo 
spin gap due to superconducting transition j20[ |2l} |22] |. 

Magnetic neutron scattering intensity is a convolution 
of instrument resolution function with dynamic magnetic 
structure factor S(ui,q). Since q scans are resolution- 
limited at Q=(7r, 7r), it is convenient to write S(uj, q) as 



S(tu,q) 



(1) 



7r 1 — exp( 

where J rfqF w (q) = ljl(|. Sharp structure of F ul (q) 
peaking at (71", w) can not be resolved in this work. Data 
in Fig. H thus represent local dynamic magnetic structure 
factor 



Sq(w) = J dqS(^,q) 



1 



7r 1 — exp(— Tuo/hsT) 



XqM, (2) 



where x'q(w) is the imaginary part of the local dynamic 
magnetic susceptibility. Relaxational spin dynamics at 
low energy is generally described byp3[ 



XqH 



(huj) 2 + r 2 



(3) 



where T is the relaxation energy and %q represents spec- 
tral intensity. For scan at each temperature in Fig. ^ T 
and xq are extracted by least squared fit of the data 
to Eqs (||) and (||), plus the usual elastic/incoherent 
peak, convoluted with instrument resolution. A flat back- 
ground measured at q=(1.39,0,0) (refer to triangle) is 
included in the fitting. The consistency between the the- 
oretical curves and measured data in Fig. ^| is very satis- 
factory. 
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FIG. 3: (a) Measured temperature dependence of F (squares) 
and Xq 1 (circles) for La2Cuo.94Lio.06O4 [Refer to Eq. (g)]. (b) 
Temperature dependence of quasielastic magnetic intensity, 
Sq(lo = + ). The filled squares are calculated using Eq. (a) 
from r and XQ m ( a )> an d open circles and diamonds were 
measured at low energy, Ttw = 0.2 meV, during warming and 
cooling, respectively. 



Experimental V and x_q as a function of tempera- 
ture are shown in Fig. @(a) by square and circle, re- 
spectively. Clearly, V and Xq 1 saturate at finite values 
at low temperatures. It is well known for conventional 
magnetic materials that in paramagnetic phase, T is a 
monotonically increasing function of T which is zero at 
magnetic transition temperature H] . This behavior ap- 
plies also to spin-freezing transition of classic spin glasses 
such as Ctt-Mnp4). When magnetic transition is pre- 
vented by frustration, the monotonic behavior for T still 
persists p5fl. Temperature independent T has been ob- 
served in intermediate- valence compounds [p6| . However, 
there is no magnetic correlations between rare-earth ions 
in these materials and the T reflects valence fluctuations. 
In La2Cuo.94Lio.06O4, magnetic correlations are substan- 
tial with £ > 42A. What we observe here is consistent 
with theoretical expectation for 2D 5=1/2 HAF with the 
QD ground statepj: T(T) ~ max[fc B T, (g - g c ) zlJ \. The 
dashed line in Fig7|(a), T(T) = max[0.18fc s T, r ] with 
To = 0.77meV, captures the main feature of our data. 

In hole-doped La2Cu04, the g is physically the dop- 
ing concentration. Besides frustrating the long-range an- 
tiferromagnetic order, doping inevitably introduces dis- 
orders. The spin freezing temperature T s f measured 
with /xSR for Li-doped La2CuC>4[^7j is nearly identical to 
that for Sr-doped La2Cu0 4 [^8|. For La2Cuo.94Lio.06O4, 
T s f = 8 K. Up to now, little has been predicted for 
S(u>,q) for realistic spin glasses. In ferromagnetic reen- 
trant spin glasses, spin dynamics expected of a ferro- 
magnet has been experimentally observed within a fi- 
nite range of q of magnetic zone center in the temper- 



ature interval T sf < T < IcQ. Along this line, one 
would expect the QD behavior to hold between ~10 and 
~60 K for La2Cuo.94Lio.06O4, refer to Fig. ||(a). In ad- 
dition, below spin freezing, T from spin clusters at (ir, n) 
in La2Cuo.94Lio.06O4 remains finite. There is no sign of 
freezing r — > 0, which is a universal feature for classic 
spin glasses. This suggests that inelastic magnetic in- 
tensity in Fig. |^ comes from spin clusters which do not 
participate in spin freezing. The upturn of r in Fig. |(a) 
below 11 K may be attributed to modification to the QD 
dynamics due to coupling of the fluctuating spin clusters 
to freezing spin clusters. 

One motivation of this low-energy neutron scattering 
work is to make connection to Cu nuclear resonance 
NMR/NQR results. Soon after the discovery of cuprate 
superconductivity, anomalous suppression of low energy 
spin fluctuations was discovered in normal state of var- 
ious family of underdoped cuprates^, |3^]. Refer to 
Ref. pi for a historic review of this so-called "pseudo 
spin-gap" phenomenon (PSG). It has been discussed var- 
iously as due to, e.g., a temperature dependent magnetic 
correlation length fl32|, |33| ], resonating valence bond sin- 
glet pairing Q, activation gaps of stripe glasses |3E|, [36[ , 
or incipient superconducting order^l). In the La2Cu04 
family, extensive inelastic neutron scattering studies have 
found no such gapp^, |l7j |, except the spin gap due to su- 
perconducting transition ]22|. 

The nuclear spin-lattice relaxation rate 1/T± is a 
weighted summation of dynamic magnetic structure fac- 
tor 5(o>,q) over the Brillouin zone at small to = + , 

Tf 1 = Kq)| 2 S(0 + , q), (4) 
q 

where hyperfine coupling |A(q)| 2 peaks at (tt, 7t) for Cu 
NMR/NQR§|. Using Eqs. (|)-(g) for 5(w,q), 

IT 1 - 5 Q (0+)^|A(q)| 2 F w (q) 
q 

~ Sq(0) = ^f. (5) 

In Fig. |(b), 5 Q (0) calculated using Eq. (fif) from the ex- 
perimentally determined xq and T in Fig. ]9(a) are shown 
together with Sq(u>) measured at Tiuj = 0.2 meV, which 
is the lowest energy transfer without significant elas- 
tic contamination (refer to Fig. ||). Circles (diamonds) 
were measured during warming (cooling) cycle. These 
data, Sq(uj ~ 0) ~ l/?i, bear remarkable resemblance 
to l/T\ observed in Sr-doped La2Cu04 of similar hole 
concentrationpl[ . The reduction of 1/T\ at low T below 
the extrapolation from high T behavior is the experimen- 
tal observation termed "pseudo spin gap" in NMR/NQR 
studies on cuprates. We have shown that there is no gap 
in spin excitation spectra for La2Cuo.94Lio.06O4. The ap- 
parent PSG behavior in Fig. |^(b) is merely a consequence 
of a saturating T at low T, which is expected for 2D 
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5=1/2 HAF of g > g c with the QD ground state. Thus a 
possible link between the NMR/NQR PSG behavior with 
the QD state is suggested by this study; specifically, the 
PSG phenomenon may not require a pseudo gap for its 
explanation. 

From this point of view, if we take doped cuprates 
without the Neel order as the experimental realization of 
g > g c , then the ubiquitous NMR/NQR PSG behavior is 
expected as a consequence of a finite r ~ {g—g c ) zv '• The 
absence of spin or pseudo spin gap from inelastic neutron 
scattering studies can then also be reconciled with the 
PSG from NMR/NQR studies. Similar NMR/NQR spec- 
tra for underdoped La2Cu04|2^, |35|, |3(| |37j] may also be 
understood as a consequence that the dopant elements, 
dopant location (in or out-of-plane) and crystal struc- 
ture (orthorhombic or tetragonal) are irrelevant pertur- 
bations to the QD regime. As such, study of the PSG 
phenomenon should pay at least as much attention to 
saturating T as to the presence of a real spin gap. 

In summary, we found dynamic spin clusters in hole- 
doped La2Cuo.94Lio.06O4 to remain (w, 7r)-correlated. 
Spin dynamics of such clusters, with correlation length 
much larger than mean impurity distance below 150 K, 
is relaxational at (jr, n). Confirming predictions of quan- 
tum critical theory for 2D 5 = 1/2 Heisenberg antifer- 
romagnet in the quantum disordered regime, the energy 
scale r of La2Cuo.94Lio.06O4 saturates around 0.8 meV 
below ~60 K. We discuss the possible connection between 
the so-called "pseudo spin-gap" phenomenon discovered 
in NMR/NQR studies on underdoped cuprates and the 
"quantum disordered" behavior observed here. In other 
words, the NMR/NQR "pseudo spin-gap" phenomenon 
may be accounted for without pseudo spin gap in spin 
excitation spectra. 
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